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Abstract. In this paper we have performed second order normalization in the generalized photogravitaional 
restricted three body problem with Poynting- Robertson drag. We have performed Birkhoff 's normalization of 
the Hamiltonian. For this we have utilised Henrard's method and expanded the coordinates of the third body 
in Double d'Alembert series. We have found the values of first and second order components. The second 
order components are obtained as solutions of the two partial differential equations. We have employed 
the first condition of KAM theorem in solving these equations. The first and second order components are 
afli'ected by radiation pressure, oblateness and P-R drag. Finaly we obtained the third order part H3 of the 
Hamiltonian in /J'^^/j''^ zero. 
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1. Introduction 

The restricted three body problem describes the motion of an infinitesimal mass moving under the gravi- 
tational effect of the two finite masses, called primaries, which move in circular orbits around their centre of 
mass on account of their mutual attraction and the infinitesimal mass not influencing the motion of the pri- 
maries. The classical restricted three body problem is generalized to include the force of radiation pressure, 
the Poynting- Robertson effect and oblateness effect. 

J. H. Poynting (1903) considered the effect of the absorption and subsequent re-emission of sunlight by 
small isolated particles in the solar system. His work was later modified by H. P. Robertson (1937) who used 
a precise relativistic treatments of the first order in the ratio of he velocity of the particle to that of light. 

The effect of radiation pressure and P-R. drag in the restricted three body problem has been studied by 
Colombo et al. (1966), Chernikov Yu. A. (1970) and Schuerman (1980) who discussed the position as well 
as the stability of the Lagrangian equilibrium points when radiation pressure, P-R drag force are included. 
Murray C. D. (1994) systematically discussed the dynamical effect of general drag in the planar circular 
restricted three body problem, Liou J. C. et al. (1995) examined the effect of radiation pressure, P-R drag 
and solar wind drag in the restricted three body problem. 

Moser's conditions (1962), Arnold's theorem (1961) and Liapunov's theorem (1956) played a significant 
role in deciding the nonlinear stability of an equilibrium point. Applying Arnold's theorem (1961), Leontovic 
(1962) examined the nonlinear stability of triangular points. Moser gave some modifications in Arnold's 
theorem. Then Deprit and Deprit (1967) investigated the nonlinear stability of triangular points by applying 
Moser's modified version of Arnold's theorem (1961). 

Bhatnagar and Hallan (1983) studied the effect of perturbations on the nonlinear stability of triangular 
points. Maciejewski and Gozdziewski (1991) described the normalization algorithms of Hamiltonian near 
an equilibrium point. Niedzielska (1994) investigated the nonlinear stability of the libration points in the 
photogravitational restricted three body problem. Mishra P. and Ishwar B.(1995) studied second order 
normalization in the generalized restricted problem of three bodies, smaller primary being an oblate spheroid. 
Ishwar B.(1997) studied nonlinear stability in the generalized restricted three body problem. 



In this paper wc have performed BirkhofF's normalization of the Hamihonian. For this we have utihsed 
Henrard's method and expanded the coordinates of the third body in Double d'Alembert series. We have 
found the values of first and second order components. The second order components are obtained as 
solutions of the two partial differential equations. We have employed the first condition of KAM theorem in 
solving these equations. The first and second order components are affected by radiation pressure, oblateness 

1/2 1/2 

and P-R drag. Finaly we obtained the third order part H3 of the Hamiltonian in I^' I2 zero. 



2. Location of Triangular Equilibrium Points 
Equations of motions are 

dUi WiNi 



(1) X — 2ny = Ux, where, Ux 



dx rf 



(2) y + 2nx = Uy, Uy = — - 

(3) ^^^nV+y^)^(l-M)gi^M^M^ 

2 ri r2 2r^ 

= (x + Ai)' + y^ ri = + l)' + y2, 71^ = 1 + ^^2, 

_ ^ ^ _ ^^"^2^^ < mi,m2 be the masses of the primaries, A2 = ^^Ji" be the oblateness 

coefEcicnt, andrp be the equatorial and polar radii respectively r be the distance between primaries, 
q = (1— be the mass reduction factor expressed in terms of the particle's radius a, density p and radiation 

pressure efficiency factor x (in the C.G.S. system) i.e., q = 1 — ^'^^ap — ~' Assumption q = constant is 
equivalent to neglecting fluctuation in the beam of solar radiation and the effect of solar radiation, the effect of 
the planet's shadow, obviously q <1. Triangular equilibrium points are given hy Ux = 0,Uy = 0, z = 0,y ^ 0, 
then we have 

r nW^i[(l-M)(l + i^2)+Ml-#)f] S^A^] 

(4) x^=xo<l 



(5) y*=yo< 



3/i(l - iJ,)yoXo 2 xo j 

nW^S^ [2m - 1 - Ml - + 7(1 - ^'(1 - f ■ 
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3/i(l-/x)?/^ y^ j 

where xq = ^ — n, yo = iS(^l — ^)^^^ and S = ql^^, as in preprint Kushvah & Ishwar (2006) 



3. Second Order Normalization 
The Lagrangian function of the problem can be written as 

(6) ^^2^^ +y )+ n{xy - xy) + —{x +y) + V2 2^ 

, i{x + ^i)x + yy y \ 

+ Wi< —^5 n arctan - ' 



2rl {x + ii). 



and the Hamiltonian \s H = —L + PxX + Pyy, where Px, Py are the momenta coordinates given by 

^ dL Wi, , ^ dL Wi 

= ai=''""^+2^^'' + ^^' a^ = ^ + "^+2;f^ 



For simplicity wc suppose qi = 1 — e, with |e| << 1 then coordinates of triangular equilibrium points can be 
written in the form 



(7) 



(8) 



- 2 _ 1 _ ^ 4. ^ 

^~ 2 3 2 +3 



6^/3 



y = 



\/3 f _ ^ _ ^ _ 2^26 ^ (1+7) 



2 I 



9V3 



27\/3 

nWi 



27\/3 



nWi 



where 7 = 1 — 2/x. We shift the origin to L4. For that, we change x 
a = + H,b = so that 



a;* + a; and y ^ y* + y- Let 



(9) 



(10) 



1 J 2e 2A2e (9 + 7) 876 

a=-<-— -A2H ^ r=r^nWi ^ 

2| 3 3 3\/3 27\/3 



\/3f 2e A2 2^26 , (I+7) 47e 



Expanding L in power series of xandy, we get 

(11) L = Lo + Li+L2 + L3 + --- 

(12) H = Ho + Hi+H2 + H3 + --- = -L + P^x + Pyy 

where Lq, Li, L2, L3 . . . arc constants, first order term, second order term, . . . respectively. Third order term 
H3 of Hamiltonian can be written as 



(13) 

where 

(14) 

(15) 
(16) 



(17) T4 



(18) 



H3 = -L3 = -- {x^Ti + 3x''yT2 + Sxy^n + y^T^ + 6T5} 



Ti = — 



16 , , 979 , 
-e + ^A2-—A2e. 



(143 + 97)^^^^(459 + 3767)^^^^ 



6V3 



27\/3 



+7 

3\/3 
1^ 



4e , 1507 , (215 + 297) 2(1174+1697) 
<^ 14 + — + 25^2 - -77^^26 - _ '' nWi - — -= ^nWie 



18 



6\/3 



27\/3 



■ , 16 A2 367 , 115(1 + 7) (959-1367) 



r32e 382 , (511 + 537) (2519- 247) „^ 

+7 <^ — + 4OA2 - —A2e + ^-^nWx - ^-^nW^e 

1 3 9 6\/3 27\/3 



1:3 = — 



9 
16 



8 203^2 625 , (105+157) „^ (403-1147) 

-e H r:r^2e 1= nWi ^ 

3 6 54 l8\/3 81V3 



nWie 



+7 



-9V3 



16 

+7 



55>l2 797 , 197 + 237 211-327 „^ 

+ — — ^2e + ^r-^nWi - ^ ^r-^nWie 

2 54 18^3 8IV3 

23^2 (37 + 7) (219 + 2537) 
— 44^26 - -^^ ^nWi - ^^-^nWie 



18\/3 



81 VS 



r. 88 

< 4eH 

\ 27 



(241 + 457) „r (1558- 1267) 
A2e + ^=—^nWi - ^ — -nWie 



18\/3 



81 V3 



n = 



2(a2 + 62)3 ~ ~ ~ + + + 

In order to perform Birkhoff's normalization, we use Henrard's method[Deprit and Deprit Brtholome 
(1967)] for which the coordinates {x,y) of infinitesimal body, to be expanded in double d'Alembert series 
x = X]n>i -^n °' y — Sn>i where the homogeneous components and B'^^ of degree n are of the 
form 



(19) 



0<m<n 



^2 X/ 



s,p,9 cos (p(j!)i + q<j)2) + sin (p(?ii + g'(j!)2) 



The condition in double summation arc (i) p runs over those integers in the interval < p < n — m, that 
have the same parity as n — m (ii) q runs over those integers in the interval —m< q < m that have the 
same parity as m. Here 7i, I2 are the action momenta coordinates which are to be taken as constants 
of integer, (f>i, (f>2 arc angle coordinates to be determined as linear functions of time in such a way that 
01 = wi + J2n>i hn{h,l2), 02 = W2 + Z]„>i g2n{h,l2) where 0Ji,uj2 are the basic frequencies, /2„ and g2n 
are of the form 



(20) 
(21) 



E£/ jn—m jm 

J 2(n-m),2m-'l ^2 



f2n — J 2(n-m),2m-'l 

0<m<n 
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EJ jn—m jm 

9 2(n-m),2m^l ^2 



0<m<n 



The first order components B\''^ and -B^'^ in 7i, I2 are the values of x and y given by 



(22) 



X = JT where X = 



X 

y 

Px 

Py 



> J — ['^u]l<»<J<4) T = 



Qi 
Pi 

P2 



(23) 



{'ZhuJi)^/'^ cos <pi, 



Q^ 



sm( 



1,2) 



(24) 
(25) 

where 

(26) 



1,0 



Jl3 ^2u!ili COS 01 + Ji4 ^2u!2h cos 02 



j2j j2J 
B?'^ = J21 -i / — - sin 01 + J22\ — - sin 02 + J23 \^^J^coi cos 0i + J24 \/2I^c02 sin 02 

V Wl V W2 



Jl3 = 



2a;ifci 



1 - 
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45A2 717^26 (67 + 197) 



36 



12\/3 



nWi 



(431 - 37) 
27\/3 



nWie 
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J_ 
2/f 

1 
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4A;2 



3^ _ 29^ _ (187 + 277)^^^ _ 



36 



3A2 



12V3 
73^26 (1 - 97) 



e - 3A2 - 



24 
299^26 
72 



+ 



24^/3 

(6 - 57) 
12v^ 



27a/3 
(53 - 397) 



+ 



2u2 



4qk 



3^2 (33 + 147) ■ 
4 I2V3 



+ 



54^/3 

(266-937) w ' 

nWi - ;=-'-^nWi€ 

54^3 

347^2 (43 - 87) 
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2j,2 



36 



4\/3 



nWi 



(27) 



Jl4 = 



2(^2^2 



2/2 



2/^2 



2/2 



717^26 



■ 293^2 , 
3e h 



45^2 

2 36 
(187 + 277) 
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-nWi p — nWie 



I2V3 



27^3 



36 



I2V3 



2(247 + 37) „/ ' 

nWi ^ ^ — -nWie 

27 v 3 



3^2 



73^26 , (1 - 97) (53 - 397) ■ 



+ 



2^2 



4/2A;2 



24 24v^ 
299^26 _ (6 - 57) 
^72 12V3 
33A2 (1643 - 937) 



e - 3^2 - 



54v^ 



216\/3 



-nWi 



7e 

4/2fc2 



54\/3 
737^2 (13 + 27) 



72 



^/3 



nWi 



Aruvi f 1 



J_ 
1 

7 



e + 



45^2 717^26 
36~ 



(67 + 19,) (413-3^^^^^- 



12V3 



27^3 



3^ _ 293^ ^ (187 + 27,)^^^ _ 2(247 + 3,)^^^^- 



36 ■ 12\/3 

73^26 (1-97) 
24 24^/3 



3^2 



nWi 



27V3 
(53 - 397) 
54^3 



nW\e 



e - 3^. - ^ - ^nW^. - i^^^nli^.e 



4n(x.'2 

2^ 

2Zi 



33^2 



72 l2\/3 

)8 - IO7) 



24v/3 



7e 



54\/3 

242^2 (43 - 87) 
9 ^ 4^3 



1 



1 

. 2l| _ 

" 293A; 
3e — — 



45^2 



717^26 

2 36 
^ (187 + 277) 
36 



+ 



1 
2^2 

7 



+ 



I2V3 

"e 73^26 (1-97) 

- - 3A2 -f- + jJ^nWi 

.2 24 24\/3 

299A2e (6-57) ^,,, 

.^-'^^ 72 WT""^^ 

■33A2 



^yyA2e ID - C7j 

WT""^^ ■ 

(34 + 57) ^^1^^ 

4Zifc| [ 4 ^ l2^/3 J 8/1^1 



I2V3 27\/3 
27V3 

(53-397) „r 

54a/3 
54v^ 

.2 4x/3 J J 



V3 r„ ^ , 37^26 
' 2e + 6vl2 + 

33^26 



4ci;iZifci 



(13 + 7) 2(79-77) TT. 



2e 

6 + y + 13^2 



2\/3 

(11 -7) 
2^/3 



9\/3 

(186 -7) 
9V3 



nWie 



+ 



1 



2Zf 

JL 

21? 

7 

7e 
'8ZfA;2 



51^2 + ^^i±SllnW^, 



3\/3 



3^^^(19 + 67)^^^ 



^ ^ 808^26 (67+197) 
6e + 135^2 7-=- - i=-^nWi 



3e - 18^2 - 



9 

55^26 



2^3 
(1 - 97) 



_ (755 + 197) 
9V^ 



nWie 



9A2 
2 



4 

(34 - 57) 
2VS 



4^3 



923-6O7 
I2V3 



\/3 37^26 
' 2e + 6^2 + 

33^26 



4:U)2l2k2 



2e 

6 + y + 13^2 - 



(13 + 7) „^ 2(79-77) 
2V3 qVs 

(11-7) (186 - 7) ■ 



1 

"2Z| 



,,^^^ (14 + 87) ^^^' 
3\/3 



e 



2/i 

7 

"2fc2 

7e 



,0.. 808^26 (67+197) (755 + 197) 
6e + I35A2 - — ^^y^nWi - ^ \. " nWie 



3e - 18^2 - 



9 

55^26 



2^3 
(1 - 97) 



9v^ 



99^2 



4 

(34 - 57) 
2V3 



4\/3 



923-6O7 „^ 
12^/3 



with Ij = 4llJj + 9,{j = 1,2) and ~ — = — 2ijj| + 1. In order to findout the second order 
components ^B^^"^ we consider Lagrange's equati 

(32) #(&-|^ = 0, 



d ^dL^ dL 
dt dy dy 







(33) 



I.e. 



X - 2ny + {2E - ■n?)x + Gy 
X + 2nx + (2F - n'^)y + Gx 



dx 



dx 



Since x and y are double d'Alembert series, x^x'^lj > 0, A; > 0, j + fc > 0) is also a double d'Alcmbert series, 
the time derivatives x,y,x,y are also double d'Alembert series. We can write 

X = ^Xn, y = ^yn, X=^X„, y = ^ j/„ 

n>l n>l n>l n>l 

where y, x, y are homogeneous components of degree n in I^^, /2^^ i-^- 



n>l n>l 

We write three components ii, X2, is of i 



9^1,0 5^1,0 

— -^(^1 + /2 + /4 + • • • ) + ^rf^(-^2 + 52 + .94 + • • • ) 



(35) 
(36) 



(37) 

where 

(38) 



Xi 



X2 



X3 



UJl- 



1,0 



,1,0 



^2^r, — = -0-°i 



dB\' 



d(^2 
dB^ 



dB 



UJ2- 



1,0 



dB 



1,0 

2 ^DBl^' 



1,0 



UJl- 



= DB'^" + f2 



LU2- 



dB{ 



1,0 



1,0 



52- 



OBl 



1,0 



D = uji- 



92- 



UJ2 



dBl 



1,0 



d(j)2 



Similarly three components Xi,X2,X3 of x are 



n2Dl O •• 7-i2 d1,0 

a;2 = B> B2 , a;3 = B> B^ 



2c^i/; 



In similar manner we can write the components of y, y. Putting the values of x, y, i, y, x and ij in terms of 
double d'Alembert series in equation f I33|l we get 

3 



(39) 



2E-1 



-Ao B 



2 1 



A2]D-G'>B 



X2 



(40) 
where 



2{1 + -A2] D + G}B2 



1,0 



X2 = 



dx 



2F - 1 - -^2 1 = 



and Y2 



9y 



These are two simultaneous partial differential equations in B^''^ and i?2'^. We solve these equations to find 
the values of B^ " and B^ ^, from Eq. ( EHl and ( EOJ 



(41) 



AiAiB2'° = $2, A1A1S2'' = -*2 where Ai = + tj^, A2 = 



nO,l 



(42) $2 = {D^ + 2F~ n^)X2 + (2ni? - G)Y2 
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(43) *2 = (2ni:> + G)X2 - (D^ + 2E- n^)Y2 

The Ea.f l41|l can be solved for B^'^ and iJj'^ by putting the formula 

cos{p(l)i + q4>2) -| f cos(p0i + 

or 



^^^"^ [ sin{p(j)i + q(j)2) '^I'^z y sin(p0i+g02) 
where 

Ap,g [ujI - [uip - ^2?)^] [i^l - (uJiP - W2q)^] 

provided A^^^ ^ 0. Since Ai^q — 0, Aq.i — the terms cos 0i, sin^i, cos 02, sin(/i2 are the critical terms. 
and tfj are free from such terms. By condition(l) of Moser's theorem fci^i + k2Uj2 7^ for all pairs (fci, fc2) 
of integers such that |fci| + |fc2| < 4, therefore each of uji,uj2,^i ± 2a;2,tf2 ± 2ti;i is different from zero 
and consequently none of the divisors Ao,o, Ao,2, A2,0! Ai^i, Ai__i is zero. The second order components 
B2' ,-82' are as follows: 



(44) bI'° = n/i +r2/2 +7-3/1 cos 2(/)i +r4/2 cos 202 +r5/y^/2^^cos(</) 

+rell^^l2^^ cos(0i + 02) + rjli sin 20i + rsh sin 202 
+rgll^'^ I^^'^ sin(0i - 02) + rio/y^/2^^ sin(0i + 02) 

and 

(45) S"'^ = - |si/i + S2/2 + S3/1 cos 201 + S4/2 cos 202 + 55/1^^/2^^ cos( 

+56/^^/2^^ cos(0i + 02) + S7/1 sin 201 + ssh sin 202 
+ sg/y^/j^^ sin(0i - 02) + Sio/y^/2^^ sin(0i + 02)| 

where 



(46) n 



f^jl^iJlFi + Ji3J23WlF^ + + J2'3^l) i^i'l 



(47) r2 = i JliU02Fi + Jl4 J24t^2Fi + ( — + J|4C^2 ) Fi" 



(48) r3 



8Ljfj2i{Ji3Fi + 2J23F[') + Au;l 



{J13F2 + J2?,F2)Ji3UJi 
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J23CJ1 F" 



-2u;iJ2i(Ji3F3 + 2J23J^3') - ^iJMsFi + J23Fi')u;i + ( ^ ~ ^23^1 j^i" 



(49) r4 



3l^|(4^2-^?) 



8cj| ,722(^14^1 + 2J24/^n - 4cc'^ 



(J14F2 + J24-F2 )Jl4:^2 ~ J24'-^2 ]F2 

W2 



-2a;2>/22(>/l4/3 + '2J2iF^') — UJ2Jia{JiaF4, + J24Fl')uj2 
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-22. — j|4a;2 j-F" 

W2 / 



(50) r5 = 



1 



(Wi + LV2f 



-2{J2lJ24(^)^/'-J22J23(-)^/'}Fi" 



{2{ Jl3 J14-F2 



+ (Jl3./24 + Jl4J23)i^2}Kw2)'/' + { 



"^'^■^'^ +J23J24(a;iU;2)l/^}F^' 



{ Jl3J22(^^)^/^ - Jl4J2l(^)^/^}F^ - 2{ J21 J24(^)^/^ - 722723(^)^/^^3" 
U>2 Wi Wi U)2 



(a;ia;2)V2 



-(wi +a;2) 

{2{ J13J14-F4 + (J13J24 + Ji4J23)i^i}('^i'^2)'/' + 2{ + ^23J24(a;ia;2)'/'}i^i' 



(51) r6 = 



(wi - ^2)^ 



t<;iCJ2(2ti-'i — a;2)(4ciJi — 2lu2) 

+ 2{ J21 J24(^)^/^ + J22 J23(^)^/^}i^l" 



{Jl3J22(^)^/^-Jl4J2l(^)^/ni^/ 
W2 ^1 



UJ2 



+(wi - a;2)^ 
-2{ 



{2{Jl3Jl4F2 + (J13J24 + Jl4J23)i^2}(wi'^2)'/' 



■iuJlOJ2y/^ 



J23J24{C01U,2)^/^}F^' 



{OJI - UJ2) 



{Jnj22{^y/' 



UJ2 



-Jl4J2l{-r^'}F^ + 2{ J2I J22(-)^/^ + J22 J23(-)^/'}F^' 
UJi U>1 U>2 



{2{ JisJw^k + (J13J24 + Ji4J23)i^i}(wia;2)'/^ - 2{ 



J21J22 
(c^ia;2)i/2 



- J23J24(a;iW2)'/'}Fi' 



(52) 



ry = 



3wf(4wf - w^) 



jso;? j^3(j^3i^^ + J23i^i')'^i - - "^2^1) 



-2a;i 



^iJlsiJlsFi + ■^23-p3) 



/2 



-4a;2j2i(Ji3i^2 + J23F^')lOi + J2l(Jl3^1 + 2J23i^i') 



(53) 



+'i0Jp22{Jl4F2 + 2J24F^')l02 " 2^2 



Jl4{Jl4Fi + J24Fl)u)2 — ( — ^ — -/24'^2 ) -F" 

W2 



-^22 



■^22 
U)2 



- ^24'^2 -F3' 



W2 Jl4('/l4-F3 + •^24-F3) 
— J22{Jl4Fl + 2J24-F4') > 



(54) 



r-9 = 



a;ia;2(2a;i + W2)(wi + 2^2) 



{uJi+UJ2f 



{2Ji3 Ji4i^i 
J21J22 



+ (Jl3J24 + Jl4J23)i^0('^1^2)'/' + / Tll2 + ■^23^24(^1^2)'/' jF^' 

(a;iCiJ2 

-(Wl + U)2f' 



{Jl3J22(^^)'/'-Jl4J2l(^^)'/'}f^ 



a;2 



-2{J21 J24(^)'/' - J22 J23(^)'/n^2' 

Wl OJ2 



-{(Jl +UJ2) 



{2{ J13J14-F3 + (J13J24 + Jl4J23)-F3}(^l'^2) 



1/2 



+2{ /'"tl/2 +-^23J24(u;iU;2)'/'}d - [{Jl3J22(^^)'/' 

^ (a;iW2) ' J L 



CU2 



- Jl4J2l(^)'/'}Fi - 2{ J21 J24(^)'/' - J22 J23(^)'/'}i^i' 



no = 



u)\UJ2{'2.u!\ — a;2)(2a;2 — Wi) 



(Wi - ljJ2f 



{2Ji3 J14F1 



(55) +(Jl3J24 + Ji4J23)i=^i'}(wia;2)'/' - 2{ 



J21J22 

(wia;2)i/2 



J23J24(a;ia;2)'/'}Fi" 



-(wi - a;2)' 

-(wi - UJ2) 



{ Jl3 J22(^)'/' - Jl4 J2l(^)'/'}i^^ + 2{ J21 J24(^^)'/' + J22 J23(^)'/'}F^' 
UJ2 ^1 ^1 ^2 



{2{ Ji3Ji4i^3 + (J13J24 + Jl4J23)F^}(u;iW2)'/' " 2{- 



J21 J22 



+ 



(a;iu;2)V2 

{ Jl3 J22(^)'/' - Jl4 J2l(^^)'/'}Fi + 2{ J21 J24(^^)'/' - J22 J23(^)'/'}F. 



J23J24(W1W2)'/'}F3' 



W2 



Wl 



Wl 



W2 



We can write expressions of s, with the help of replacing Fj by G,, F- by G- and i^/' by G'!,{i = 1, 2, 3, 4) 
where 



(56) 
(57) 



(59) 



F2 



-nWie 



3 

32 



16 ^ , 979 , (143 + 97) (555 + 3767) 



4e 



+7 <; 14 + y + 25^2 



1507 , (215 + 297) „^ 2(1174+1697) 



(58) F3 = 



3V3 
1^ 



, , 16 23^2 

14 eH ^ 

3 2 



18 

104 



A2e + 



6^/3 
115(1 + 7) 



27V3 
_ 2(439- 687) 
27\/3 



r32e ^„ , 310, (511 + 537) „^ (2519- 2497 

+7 <^ ^ + 40^2 - ^^2e + ^-^nWi - ^^nW'ie 

I 3 9 6^/3 27V3 



-3 
256 



o^. 17801^2 , (2821 + 1897) (23077 + 95927) „^ 
364 + 420^2 7: — -A^e + -= '-^nWi - ^ —= '-!-nWie 



3\/3 



27\/3 



r ,^ , lOOe , 849^2 , 59 ^ ^ (125 + 387)^^^^ (87613- 2137) 



+287|23+ 2^ 



14 



+ y ^26 • 



6\/3 



27\/3 



(60) Fi 
(61) 



(63) Fi 



3^3 
3^3 
1^ 



14- 



16 



1367 , 115(1 + 7) (863- 1367) 



18 



isVs 



r32e , 382 , (511 + 537) 



27\/3 

(2519 - 247) 
27V3 



nW-] e 



(62) = — 



8 203^2 721 (105 + 167) (319-1147) 

-e H ;7-r^2e - -^^ 1= — -nWi - ^ — -nWie 

3 6 54 ^ 18V3 8IV3 



2-1! 
9 



173Ao 



781 , 

9 ^ 



(197 + 237) 



-3V3 
16 



392 



6 9 "^^ ' 18V3 

532e 1918^2 28582^2 



(265 - 327) 
81V3 



+287 <! + 



3 3 

4037^2 611 



, , (203 + 12II7) , (949 + 43787) 
A2e + ^ F ^nW^i + -^^ 7= — -nWie 



9^3 



27V3 



84 



(8397 + 9197) (92266- 18697) 

-nVVi 7= nlTie 



84\/3 



27\/3 



(64) Fi' = 

(65) F^' = 

(66) F^' = 

(67) Fi' = 



nWic 
6 

-9 " 
32 



+7 



— e 
3 



203^2 
6 

4e 55A2 



625 , (105 + 157) (307-1147) 
54 18V3 8lV3 

797 , (197 + 237) (211-327) „^ 

54 18V3 8IV3 



55^2 
6 



134 



-A2e 



(37 + 7) 
18^3 



(93 + 2267) 
81V3 



-9V3 
16 

( 169, (241+457) TT. (1558- 1267) „^ 

\ 27 18x/3 81\/3 



nWie 



9 
256 



212 295OA2 1370^2 , (771 + 2377) „^ 

-A2€ ■= nWi 



27 



9V3 



2(1907 - 9847) 
8IV3 



nWie 



36965 , (2569 + 2777) (22603 + 43967) 

-A2e + -j^-^nWi + ■j= — -nWie 



504 



252V3 



1134 V3 



(68) 
(69) 

(70) 

(71) 



Gi = 

G2 = 

G3 = 

G4 = 



-nWie 



_3_ 

32 



3V3 



14 



16 



Ao 



1367 
1^ 



^26- 



115(1 + 7) 
18V3 



nWi 



(863 - 1367) 
27^/3 



nWie 



16 

+7 



32e 382 , (511 + 537) (2519-247) 

i?e + 6^2 - ^ A2e + iil^nl^, + liII±«nW^.e 



3\/3 
256 



18 ' 6x/3 
, 71^2 1489 , (215 + 297) 

+ ^71— 

84 + 52e + 212A2 - 267^26 + ^(^^9 + 617) 

3\/3 



27\/3 

2(1174+1697) 
27\/3 
(14854 + 2257) 



27\/3 



nWie 



+7 ^ 32e + 156A2 + 649A2e - ^^^^nW, + (13285 + 51697)^^^^ 



3\/3 
10 



27\/3 



(72) 
(73) 



g; = 4 



203^2 625 
6 M 
55^2 



(105 + 157) 



nWi 



(307-1147) 



797 , 

54 ^ 



18\/3 
(197 + 237) 



(74) G'^ 



(75) 



3\/3 



1. 16 



65^2 1439 



+7 

-9 
128 



12e - 287^2 
152e 



6 

310 , 

^26 

9 ^ 



18V3 
115(1 + 7) 



18 ' 18V3 
(511 + 537) 



nWi - 
nWi 



81\/3 

(211 - 327) 
81V3 

(941 - II87) 



(76) 
(77) 

(78) 
(79) 



7<^ 96 + 



-nWie 



6 

9V3 
"32" 



847v42 
135^2 - 

23^2 



A2e- 

2320 
~9~ 



6^3 
2(28 + 7) 



nWi 



27^3 
(251 - 247) 



nWie 

nWie 



27\/3 



nW^e 



4(2210- 697) 
nWi ^ j= -nWie 



Aye 



V3 

(497 - 1237) 



27V3 

4(17697 + 327) 

nWi ^ ^ -nWiC 

27%/3 



- 44^26 



(37 + 7) (123 + 3497) 

-nWi ^ nWie 



+7 i 4e + 



88A2 (421 + 457) 



G' = ^ 



9 
16 



27 
203^2 



18a/3 



18\/3 3\/3 
(1558-1267) 



81\/3 



-nWie 



6 



589 , 5(51 + 27) 



(349 - 2827) 



g: 



2 

-9V3 
256 

+7|8e 



4e 



26^2 



412 
"27 



-Aye 



(197 + 237) 



20 

12 + -. 
749A2 



76.12 
808 



350yl, 



18^3 

(327) 



nWi 



81%/3 
(211 - 327) 
81\/3 



nWie 



nWie 



9 



Aye 



'-Aye , ^ 
3 3V3 

(IO9-4O7) 



2(1529 + 4507) ^. . 
riM/i ■= nWie 



3V3 



27V3 
(35 - 12697) 



27V3 



nWie 



4. Conclusion 



Using transformation x = bI''^ + B^''^ and y = B^'^ + Bj'^ in Eg.f I13|l the third order part of the 



rl/2rl/2 



is of the form 



ri/2 



Ai,yiy^)l2 + Ao,34^^ 



Hamihonian in Ii ly 

(80) i/3 = ^3.0/?^' + Ay^JJ^ 

We can verify that in Eq.( I80|l A3.0 vanishes independently as in Deprit and Deprit Barhtolome(1967). 

Similarly the other coefficients ^2,1 j ^1,2; ^0,3 Etre also found to be zero independently. Hence the third 

1/2 1/2 

order part i?3 of the Hamiltonian in I-^ I2 is zero. 
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